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Abstract 

This paper studies the geodesic diameter of polygonal domains having h holes and n corners. 
For simple polygons (i.e., h — 0), it is known that the geodesic diameter is determined by a pair of 
corners of a given polygon and can be computed in linear time. For general polygonal domains with 
h > 1, however, no algorithm for computing the geodesic diameter was known prior to this paper. 
In this paper, we present first algorithms that compute the geodesic diameter of a given polygonal 
domain in worst-case time 0(ri^-^'^) or 0{n7 {log n + h)). The algorithms are based on our new 
geometric observations, part of which states as follows: the geodesic diameter of a polygonal domain 
can be determined by two points in its interior, and in that case there are at least five shortest paths 
between the two points. 
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Science and Culture, Japan, and Japan Society for the Promotion of Science. 



1 Introduction 



A polygonal domain V with h holes and n corners 1/ is a connected and closed subset of of genus 
h whose boundary &P consists of /i + 1 simple closed polygonal chains of n total line segments. The 
holes and the outer boundary of V are regarded as obstacles so that any feasible path in V is not allowed 
to cross the boundary dV. The geodesic distance d(p, q) between any two points p, g in a polygonal 
domain V is defined as the (Euclidean) length of a shortest obstacle-avoiding path between -p and q. 

In this paper, we address the geodesic diameter problem in polygonal domains. The geodesic di- 
ameter diam(7') of a polygonal domain V is defined as diam(T') := max^ d(s, t). A pair (s, t) 
of points in V that realizes the geodesic diameter diam(T') is called a diametral pair. The geodesic 
diameter problem is to find the value of diam('P) and a diametral pair. 

For simple polygons (i.e., h = 0), the geodesic diameter has been extensively studied and fully 
understood. Chazelle ||6l provided the first O(n^) -time algorithm computing the geodesic diameter 
of a simple polygon, and Suri fl8| presented an 0(n log n) -time algorithm that solves the all-geodesic- 
farthest neighbors problem, computing the farthest neighbor of every corner and thus finding the geodesic 
diameter. At last, Hershberger and Suri lITTll showed that the diameter can be computed in linear time 
using their fast matrix search technique. 

On the other hand, to the best of our knowledge, no algorithm for computing diam(T') has yet been 
discovered when P is a polygonal domain having one or more holes {h > 1). Mitchell [14J has posed 
an open problem asking an algorithm for computing the geodesic diameter diam('P). However, even 
for the corner-to-corner diameter max^ d(u, v), only known is a brute-force algorithm that takes 
0(n^ log n) time, checking all the geodesic distances between every pair of cornersjj 

This fairly wide gap between simple polygons and polygonal domains is seemingly due to the 
uniqueness of the shortest path between any two points; it is well known that there is a unique shortest 
path between any two points in a simple polygon |l9l- Using this uniqueness, one can show that the diam- 
eter is indeed realized by a pair of corners in V; that is, diam('P) = max^ ,;gy d('u, v) if h = II11U18II . 
For general polygonal domains with h > 1, however, this is not the case. In this paper, we exhibit 
several examples where the diameter is realized by non-comer points on dV or even by interior points 
of V. (See Figure [Hand Appendix A.) This observation also shows an immediate difficulty in devising 
any exhaustive algorithm since the search space like dV or the whole domain V is not discrete. 

The status of the geodesic center problem is also similar. The geodesic center is defined to be a 
point in V that minimizes the maximum geodesic distance from it to any other point of V. Asano and 
Toussaint IS introduced the first 0(n^ log n)-time algorithm for computing the geodesic center of a 
simple polygon, and Pollack, Sharir and Rote fTT] improved it to 0{n log n) time. As with the diameter 
problem, there is no known algorithm for general polygonal domains. Note that computing the geodesic 
center involves computing the geodesic diameter because the geodesic center may be determined by the 
midpoint of a shortest path defining the geodesic diameter. See O'Rourke and Suri |[T6l and Mitchell lfT4l 
for more references on the geodesic diameter/center problem in simple polygons and polygonal domains. 

Since the geodesic diameter/center of a simple polygon is determined by its corners, one can exploit 
the geodesic farthest-site Voronoi diagram of the comers V to compute the diameter/center, which can 
be built in O(nlogn) time ||2l. Recently, Bae and Chwa Q presented an 0(nA;log^(n + A;))-time 
algorithm for computing the geodesic farthest-site Voronoi diagram of k sites in a general polygonal 
domain. This can be used to compute the geodesic diameter maxp d{p, q) of & finite set S of points 
in V, but cannot be exploited for computing diam(P) without any characterization of the geodesic 
diameter of polygonal domains with h>l. Moreover, when S = V, this approach is no better than the 
brute-force 0(n^ logn)-time algorithm for computing the corner-to-corner diameter m.ayiu,v&v d(u, v). 

In this paper, we present the first algorithms that compute the geodesic diameter of a given polygonal 
domain in 0(n^'^^) or 0(n^(log n + h)) time in the worst case. We also show that for small constant h 

'Personal communication with Mitchell. 
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Figure 1: Three polygonal domains where the geodesic diameter is determined by a pair (s*, t*) of non-corner 
points; Gray-shaded regions depict the interior of the holes and dark gray segments depict the boundary dV. Recall 
that V, as a set, contains its boundary &P. (a) Both s* and t* lie on dV. There are three shortest paths between 
s* and t* . In this polygonal domain, there are two (symmetric) diametral pairs, (b) s* e dV \ V and t* e intT'. 
Three triangular holes are placed in a symmetric way. There are four shortest paths between s* and t* . (c) Both s* 
and t* lie in the interior uiiV. Here, the five holes are packed like jigsaw puzzle pieces, forming narrow corridors 
(dark gray paths) and two empty, regular triangles. Observe that d(ui, wi) = d(ui, W2) = d(M2, V2) — d(u2, ^3) 
— d(u3,W3) = d(u3,wi). s* and t* lie at the centers of the triangles formed by the Ui and the Vi, respectively. 
There are six shortest paths between s* and t* . More details on this example can be found in Appendix A. 2. 

the diameter can be computed much faster. Our new geometric results underlying the algorithms indeed 
show that the existence of any diametral pair consisting of non-corner points implies multiple shortest 
paths between the pair; as one of the cases, it is shown that if (s, t) is a diametral pair and both s and t 
lie in the interior ofV, then there are at least five shortest paths between s and t. 

Some analogies between polygonal domains and convex poly topes in M'^ can be seen. O'Rourke and 
Schevon [15,1 proved that if the geodesic diameter on a convex 3-polytope is realized by two non-corner 
points, at least five shortest paths exist between the two. Based on this observation, they presented an 
0(n^'^ log n) -time algorithm for computing the geodesic diameter on a convex 3-polytope. Afterwards, 
the time bound has been improved to 0(n^ log n) by Agarwal et al. |T| and recently to 0{'n7 log n) by 
Cook IV and Wenk [8|. This is also compared with the running time of our algorithms for polygonal 
domains, 0(n''-^^) or 0(n'' (log n + h)). 

2 Preliminaries 

Throughout the paper, we frequently use several topological concepts such as open and closed subsets, 
neighborhoods, and the boundary dA and the interior iniA of a set A; unless stated otherwise, all of 
them are supposed to be derived with respect to the standard topology on with the Euclidean norm 
II • II for fixed d>l. We denote the straight line segment joining two points a, 6 by ah. 

We are given as input a polygonal domain V with h holes and n corners. More precisely, V consists 
of an outer simple polygon in the plane and a set of h (> 0) disjoint simple polygons inside the 
outer polygon. As a subset of M^, V is the region contained in its outer polygon excluding the interior 
of the holes; thus is a bounded, closed subset of M^. The boundary dV of V is regarded as a series of 
obstacles so that any feasible path inside V is not allowed to cross dV. Note that some portion or the 
whole of a feasible path may go along the boundary dV. The length of a path is the sum of the Euclidean 
lengths of its segments. It is well known from earlier work that there always exists a shortest (feasible) 
path between any two points p,q G V [ 13|. The geodesic distance, denoted by d(p, q), is then defined 
to be the length of a shortest path between p and q^V. 

Shortest path map. Let V be the set of all comers of V and 7r(s, t) be a shortest path between s £ V 
and t £ V. Then, it is represented as a sequence 7r(s, t) = (s, t^i, . . . , Vk,t) for some vi, . . . ,Vk € V; 
that is, a polygonal chain through a sequence of corners |[T3l . Note that possibly we may have A; = 
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when d(s, t) = \\s — t\\. If two paths (with possibly different endpoints) induce the same sequence of 
corners, then they are said to have the same combinatorial structure. 

The shortest path map SPM(s) for a fixed s G "P is a decomposition of V into cells such that every 
point in a common cell can be reached from s by shortest paths of the same combinatorial structure. 
Each cell as{v) of SPM(s) is associated with a corner v £ V or s itself, which is the last corner of 
7r{s,t) for any t in the cell as{v). In particular, the cell cTs(s) is the set of points t such that 7r(s,t) 
passes through no corner in V and thus d(s,t) = ||s — t\\. Each edge of SPM(s) is an arc on the 
boundary of two incident cells as{vi) and (7s{v2) and thus determined by two corners ^1,^2 G ^ U {s}. 
Similarly, each vertex of SPM(s) is determined by at least three corners vi,V2,V3 S F U {s}. Note that 
for fixed s £ V a. point t that locally maximizes ds{t) := d(s, t) lies at either (1) a vertex of SPM(s), 
(2) an intersection between the boundary dV and an edge of SPM(s), or (3) a corner in V. 

The shortest path map SPM(s) has 0{n) complexity can be computed in 0(n log n) time using 
0{n log n) working space ||T2|| . For more details on shortest path maps, see |[T2]4T4l . 



Path-length function. If tt{s, t) / st, then there are two comers u,v £ V such that tt{s, t) is formed 
as the union of a shortest path from u to f and two segments su and vt. Note that u and v are not 
necessarily distinct. In order to realize such a path, we assert that s is visible from u and t is visible 
from v; thus, s E VP(u) and t E \/P{v), where VP(p) for any p E "P is defined to be the set of all points 
q £ V such that pq C V. The set VP(p) is also called the visibility profile of p £ V [7 ]. 

We now define the path-length fi^nction leuu^v '■ VP(ii) x VP(t;) — )■ R for any fixed pair of corners 
u,v £ V tohe 

len„^^(s, t) := \\s — u\\ + d{u, v) + ||f — 

Then, leuu^vis, t) represents the length of the path from s to t that has the fixed combinatorial structure, 
entering u from s and exiting v to t. Also, unless d(s,t) = ||s — t|| (equivalently, s E VP(t)), the 
geodesic distance d(s, t) can be expressed as the pointwise minimum of some path-length functions: 

d(s,t) = mill lenuv(s,t). 
«gVP(s), ve\/P{t) 

Consequently, we have two possibilities for a diametral pair (s*,t*); either we have d(s*,t*) = 
1 1 s* — t*\\ or the pair (s* , t* ) is a local maximum of the lower envelope of several path-length functions. 



3 Local Maxima on the Lower Envelope of Convex Functions 

In this section, we give some analysis on local maxima of the lower envelope of convex functions, which 
provides a key observation for our further discussions on the geodesic diameter and diametral pairs. 

We start with a basic observation on the intersection of hemispheres on a unit sphere in the d- 
dimensional space W^. For any fixed positive integer d, let S'^~^ := {x E M*^ | ||x|| = 1} be the 
unit sphere in M'^ centered at the origin. A closed (or open) hemisphere on S'^~^ is defined to be the 
intersection of S'^~^ and a closed (open, respectively) half-space of bounded by a hyperplane that 
contains the origin. We call a fc-dimensional affine subspace of a k-flat. Note that a hyperplane in 

is a (d — l)-flat and a line in is a 1-flat. Also, the intersection of S'^~^ and a A;-flat through the 
origin in is called a great {k — l)-sphere on S'^^^. Note that a great 1-sphere is called a great circle 
and a great 0-sphere consists of two antipodal points. Then, we observe the following. 

Lemma 1 For any two positive integers d and m < d, a set of any m closed hemispheres on S'^^^ has 
a nonempty common intersection. Moreover, if the intersection has an empty interior relative to S'^^^, 
then it includes a great {d — m) -sphere on S'^~^. 

Proof. Proof can be found in Appendix B. B] 
Using Lemma [T] we prove the following theorem, which is the goal of this section. 
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Theorem 1 Let T be a finite family of real-valued convex functions defined on an open and convex 
subset C C M'^ and g{x) := minygj- /(x) be their pointwise minimum. Suppose that g attains a 
local maximum at x* G C and there are exactly m functions /i, . . . , /m G T such that m < d and 
fi{x*) = g{x*) for all i = I, . . . ,m. If none of the fi attains a local minimum at x* , then there exists a 
{d + 1 — m)-flat Lp <ZW^ through x* such that g is constant on ipHU for some neighborhood [/ C M*^ 
of X* with U C C. 

Proof. Let x* € C and m be as in the statement. For each i, consider the sublevel set Lj := {x G C | 
fiix) < fi{x*)}. Since each fi is convex and x* does not minimize fi, the set Lj is convex and x* Ues 
on the boundary dLi of Li. Therefore, there exists a supporting hyperplane hi to Li at x*. Denote by hf 
the closed half-space that is bounded by hi and does not contain Lj. Note that fi{x*) < fi{x) for any 
X G hf DC and fi{x*) < for any x G {hf\hi)nC. Let Hi := {x-x* | x G hf,\\x-x*\\ = 1} 
be a closed hemisphere on the unit sphere S'^^^ centered at the origin. 

Since g{x*) = fi{x*) for any i G {1, . . . , m} and x* is a local maximum of g, the intersection p| Hi 
has an empty interior relative to S"^^^; otherwise, there exists y G S'^~^ such that fi{x* + Ay) > ffx*) 
for any i G {!,..., m} and any A > with x* + Ay G C. Hence, by Lemma [T] p| Hi has a nonempty 
intersection including a great {d — m)-sphere G on S"^^^. Let ip be the corresponding {d — m + l)-flat 
in through x* defined as ip := {x* + Ay G M'^ | y G G and A G M}. Consider the restriction /i|(^nc 
of fi on if n C. Since fi is convex and if is an affine subspace (thus convex), fi\^nc is also convex 
and their pointwise minimum g\^nc attains a local maximum at x*. Furthermore, each fil^nc attains a 
local minimum at x*; since ip C hf, we have fi{x*) < /(x) for any point x G n C. Hence, y|(^nc 
also attains a local minimum at x* since y(x*) = fi{x*) for any i G {1, . . . , m}. Consequently, g is 
locally constant at x* on ip; more precisely, there is a sufficiently small neighborhood U C M'^ of x* 
with U C C such that g is constant on U Cup, completing the proof. 

Remark that the theorem should have its own interest and find an application in problems of maxi- 
mizing the pointwise minimum of several convex functions. 

4 Properties of Geodesic-Maximal Pairs 

We call a pair {s* ,t*) £ V xV maximal if {s*,t*) is a local maximum of the geodesic distance function 
d. That is, (s*, t*) is maximal if and only if there are two neighborhoods Us,Ut C of s* and of t*, 
respectively, such that for any s £ UsHV and any t £ UtHV we have d{s*,t*) > d{s, t). For any pair 
(s, t), let n(s, t) = {tti, . . . , TTm} bc the set of all distinct shortest paths from s to t, where m denotes 
the number of shortest paths. Let Ui and Vi be the first and the last corners in V along tTj from s to t, and 
let Vs := {ui, . . .,Um} and Vt := {vi, . . .,Vm}- 

Let E be the set of all sides of V without their endpoints and B be their union. Note that B = dV \ V, 
the boundary of V except the corners V . The goal of this section is to prove the following theorem, which 
is the main geometric result of this paper. 

Theorem 2 Suppose that {s*,t*) is a maximal pair in V andl[{s* ,t*), Vg*, and Vt* be defined as above. 
Then, we have the following implications. 



(VV) 


s* 


G V, 


t* 


G V 


implies 


\U{s* 


t*) 


>1, 


\Vs*\ 


> 1, 


\Vt* 


> 1 


(VB) 


s* 


G V, 


t* 


G B 


implies 


\U{s* 


t*) 


>2, 


\Vs*\ 


> 1, 


\Vt* 


> 2 


(VI) 


s* 


G V, 


t* 


G intV 


implies 


\U{s* 


t*) 


>3, 


\Vs*\ 


> 1, 


\Vt* 


> 3 


(BB) 


s* 


G B, 


t* 


G B 


implies 


\U{s* 


t*) 


>3, 


\Vs*\ 


>2, 


\Vt* 


> 2 


(BI) 


s* 


G B, 


t* 


G intV 


implies 


\U{s* 


t*) 


>4, 


\Vs-A 


>2, 


\Vt* 


> 3 


(II) 


s* 


G intV, t* 


G intV 


implies 


\U{s* 


t*) 


>5, 


\Vs* 


>3, 




> 3 



Moreover, each of the above bounds is best possible by examples. 
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Figure 2: (a) How to determine u'^. (left to right) Ui — u'^, s*, Ui, and the second corner are coUinear; s* and the 
first three corners are collinear (b) For points in a small disk B centered at s* with B C VP(w^) U VP(ui), the 
function /if measures the length of the shortest path from u[ to each. 

To see the tightness of the bounds, we present examples with remarks in Figure [T] and Appendix A. 
In particular, one can easily see the tightness of the bounds on \ Vs* \ and \ Vf \ from shortest path maps 
SPM(s*) and SPM(t*), when V U {s* ,t*} is in general position. 

Throughout this section, for easy discussion, we assume that there is a unique shortest path between 
any two comers u,v ^ V . This assumption does not affect Theorem |2] since multiple shortest paths 
between comers in V can only increase |n(s*,t*)|. Note that this assumption implies that the pairs 
{ui,Vi) are distinct, while the Ui (also the Vi) are not necessarily distinct. We thus have \ Vs*\ < m, 
I Vf* I < m, and \{{ui,Vi) | 1 < i < m}\ = m, where m = |n(s*, 

The following lemma proves the bounds on | Vs* \ and | Vt* \ of Theorem |2l Proofs of the lemmas 
presented in this section can be found in Appendix B. 

Lemma 2 Let (s*, t*) be a maximal pair. Then, \Vf \ > 2 ift* € B; \Vf \ > 3 ift* G iniV. Moreover, 
if t* £ e £ E, then there exists v G Vt* such that v is off the line supporting e; if t* € intT^, then t* Ues 
in the interior of the convex hull ofVt* ■ 

Lemma |2] immediately implies the lower bound on |n(s*,t*)| when s* £ V or t* £ V since 
|n(s*,r)| > max{|T4.|, \Vf\]. This finishes the proof for Cases (V-). Note that Case (VV) is trivial. 

From now on, we assume that both s* and t* are not corners in V . This assumption, together with 
Lemma[2l implies multiple shortest paths between s* and t* , and thus d{s* ,t*) > — Hence, as 
discussed in Section|2l any maximal pair falling into one of Cases (BB), (BI), and (II) appears as a local 
maximum of the lower envelope of some path-length functions. 

Case (II): When both s* and t* lie in intV. We will apply Theorem[T]to prove Theorem|2]for Case 
(II). For the purpose, we find m = |n(s*, convex functions fi defined on a convex neighborhood 
C of such that the following requirements are satisfied: (i) the pointwise minimum g of the 

fi coincides with the geodesic distance d on C, (ii) fi{s*,t*) = g{s*,t*) = d{s*,t*) for any i £ 
{!,... ,m}, (iii) g attains a local maximum at (s*,t*) G C, and (iv) none of the fi attains a local 
minimum at (s*, t*). 

If there are exactly m pairs (u, v) of corners such that leriu t*) = d(s*, t*), then we can apply 
Theorem [T] simply with the m path-length functions len^.^^,-. Unfortunately, this is not always the case; 
a single shortest path vTj G n(s*, t*) may result in several pairs (n, v) of comers with u,v £ iTi such that 
(n, v) / (uj, Vi) and leuu^vis* ,t*) = d{s*,t*). This happens only when either u, Ui, s* or v, Vi,t* are 
collinear. In this degenerate case, the path-length functions len^.^t,; violate the first requirement above. 
In the following, we thus define the merged path-length functions that satisfy all the requirements even 
under the degenerate case. 

Recall that the combinatorial structure of each shortest path vTj £ Il{s*,t*) can be represented by 
a sequence (uj = Uj i, . . . , ^ = Vi) of comers in V. We define u[ to be one of the Uij as follows: 
If s* does not lie on the line ^ C through Ui and Uj 2> then n- := n^; otherwise, if s* £ i, then 
u'i := Uij, where j is the largest index such that for any open neighborhood U C of s* there exists 
a point s £ {U n VP(nj j)) \ £. Note that such always exists, and if no three of V are collinear, then 
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Figure 3: Illustration to Lemma|4j for any {s,t) S Di and any sufficiently small S if we pick {s',t') such that 
s' is S closer to Ui than s and t' is 6 farther from Vi than t, then we have fi{s',t') ~ fi{s,t). Symmetrically, 

with lis" = \\s ~ U,\\ + S SiUd \\t" - V,\\ - \\t-V,\\-6. 

we always have either u'^ = ui or u[ = Ui^2', Figurel^a) illustrates how to determine u[. Also, we define 
v'j^ in an analogous way. Let hf and hj be two functions defined as 



where A := (VP(uO U VP(ni)) x (VP(t;0 U VP(i;i)) C P x P. We consider P x P as a subset of 
and each pair (s, G "P x "P as a point in M^. Also, we denote by {sx, Sy) the coordinates of a 

point s E "P and we write s = {sx, Sy) or (s, t) = {sx, Sy, tx,ty) by an abuse of notation. Observe that 

fi{s,t) = min{len„^,^^(s,t), len^^^^(s,i),len„^_^/(s,t),len„/^^/(s,t)} for any (s,t) G A if we define 

lenu_„(s, t) = oo when s VP(n) or v ^ \/P{v); see Figurelltb). 
We first show the convexity of the functions /j. 

Lemma 3 For any i G {1, . . . , m} and any convex subset C C Di, fi is convex on C. 

Observe that each of the /,; is indeed not strictly convex. Figure |3] illustrates one such line in for 
a fixed (s, t) G Di that fi stays constant when (s, t) moves locally along the line. We show that such a 
line in is unique for any fixed (s, i) G D^. 

Lemma 4 For any i G {1, . . . , m} and any (s, t) G intDj, there exists a unique line £i C through 
(s, t) such that fi is constant on iiHU for some neighborhood U of (s, t) with U C Di. Moreover, fi 
is constant oniiCiC for any convex neighborhood C of {s, t) with C d Di. 

Now, we let (7 : P| A — M be the pointwise minimum of the fi defined as g{s, t) = minj /j(s, t) for 
any (s, t) G H Di. Note that the intersection P| Di contains a nonempty interior and (s*, t*) G int P| Di 
by our construction. We show that the fi satisfy the aforementioned requirements to apply Theorem [T] 

Lemma 5 The functions fi and their pointwise minimum g satisfy the following conditions. 

(i) There exists a convex neighborhood C CM."^ of {s*, t*) with C C. f]Di such that d(s, t) = g{s, t) 
for any {s,t) G C. 

(ii) fiis*,t*)=g{s*,t*) = d{s*,t*)foranyi G {l,...,m}. 
(Hi) g attains a local maximum at {s*,t*). 

(iv) None of the fi attains a local minimum at {s*,t*). 

Now, we take a convex neighborhood C C of (s* , t * ) that is as described in Lemma|5] We restrict 
fii ■ ■ ■ 1 fm and g on C. Then, each fi is convex by Lemma [3] and, by Lemma |5l the m functions fi and 
their pointwise minimum g satisfy the conditions of Theorem [T] with open convex domain C C M^. 

Suppose that m < 5. Then, Theorem [T] implies that there exists at least one line i CR'^ through 
{s*,t*) such that g is constant on £r\C. On the other hand. Lemma |4] implies that there is a unique line 





Then, the merged path-length function /j : — )• M is defined as 



f,is,t) ■.= htis) + diu'„v'i) + hi{t) 



li through {s* ,t*) such that fi is constant on n C, and no other line along which fi is locally constant 
at {s* ,t*). Hence, there exists at least one index i such that I = £j. We also observe the following. 

Lemma 6 There are at most two indices i S {1, . . . , m} suctt that i = £i. 

The last task is to check two possibilities; one or two of the fi are constant on £ n C Also, recall 
that m > 3 by Lemma [2l Without loss of generality, we first assume that I = i\ ^ £i for any i > 2. 
Along l\, each fi with z > 2 is not constant but convex. Since t) = (^(s, t) = mirij fi{s, t) for any 
(s, t) S ^ n C, by LemmalU fi with i>2 must strictly increase from (s*, t*) in both directions along I. 
Thus, for any (s, t) G ^nC with (s, t) / (s*, t*), we have a strict inequality g{s, t) = t) < fi{s, t). 
Then, by Lemmas [3] and |4] at any such (s, t) G ^ n C there is a direction in which /i strictly increases: 
more precisely, for any arbitrarily small neighborhood U G C of (s*, t*), g{s, t) = t) < fi{s, t) 
for (s, t) £ dU n i and thus there exist a sufficiently small neighborhood U' C C of (s, t) and (s', t') G 
[/' such that fi{s,t) < fi{s',t') < fi{s',t') for any i > 2, which implies that g{s*,t*) = g{s,t) < 
g{s', t'), a contradiction to that g attains a local maximum at {s*,t*). 

Thus, two of the fi must be constant on ^ n C. We assume that £ = ii = £2 7^ f-i for i > 3. In 
this case, for any (s, i) £ £(1 C with (s, t) 7^ (s*, t*), we have a strict inequality ^(s, t) = /i(s, t) = 
/2(s, t) < fi{s, t). Then there exists a direction from (s, t) in which both of /i and /2 strictly increase 
by Lemmas [3] and m Since ^(s, t) = fi{s, t) = /2(s, t) < fi{s, t) for any i > 3, we get a contradiction 
analogously to the above. 

Hence, we achieve a bound m = |n(s*, > 5, as claimed in Case (II) of Theorem|2] 

Case (BB): When both s* and t* lie on B. In this case, we assume that s* G G and t* e et £ E. 
Let p be an endpoint of eg and Ig be the length of e^. We denote by s{Cs) the unique point on such that 
ll*(Cs) ~ pII = Cs for any < < h- Here, we consider s : (0, Ig) — )• as a bijective map between a 
real open interval (0, Ig) C M and a segment Cg C except its endpoints. Analogously, we also define 
t{Ct)- Let (g and Q be real numbers such that s* = s{Cs) and t* = t{Q). 

The outline of proof is analogous to the above discussion for Case (II). We redefine /j : — )• M as 

fiiCs, Ct) := KisiCg)) + d{ul v'i) + hiitiCt)), 

where A := s-i((VP(n^ U VP('Ui)) n e,) x t-^{(\/P{v'{) U VP(ui)) n e*). We consider A as a 
subset of and each pair {(^g, Q) G A as a point in M^. Also, let g{Cs, Ct) '■= min, fi{Cs, Ct) for any 

(Cs,a) eflA. 

The convexity of fi on any convex subset of Di is deduced from Lemma|3] Analogously to Lemma[5l 
one can show that (i) there exists a convex neighborhood C C of (CiO with C C f]Di such 
that5(C.,Ci) = mmfiiCgXt) = d{s{Cg),t{Ct)) for any (G,Ct) G C, (ii) fiiCCt) = diCCt) = 
d{s{Cg),t{C*)) for any i G {1, . . . , m}, (iii) g attains a local maximum at {C*, Ct)- Also, observe that 
(iv) none of the fi attains a local minimum at {Cg, Ct) if m < 3: Assume that some fi attains a local 
minimum at {C*, Ct) and m < 3. This happens only when s* = s{Cg) is the perpendicular foot of Ui on 
eg and t* = t{C*) is the perpendicular foot of Vi on et- In this case, there always exists a direction along 
eg such that if we move s* in the direction, then hf strictly increases for every 1 < i < m < 3, which 
contradicts to the assumption that {s*,t*) is maximal. 

In addition, we observe the following. 

Lemma 7 If there exists a line ^ C such that fi is constant on inC, then Ui lies on the line supporting 
eg and Vi lies on the line supporting et- 

We assume that m < 3, and restrict the functions fi and gonC C. M^. Then, Theorem[T]implies that 
there exists a line ^ C through (C^ , Ct) ^uch that g is constant on C. Moreover, Lemma |4] implies 
that at least one of the fi is constant on lr\C. Assume that only /i is constant on lr\C. Then, by the same 
argument as above for Case (II), for any point {Cs,Ct) G £r\C with (Cs, Ct) 7^ (CI) Ct*)' we have a strict 
inequality g{Cs,Ct) = fiiCsXt) < fi{Cs,Ct) for i > 2, leading to a contradiction: for any arbitrarily 
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small neighborhood U (Z C of {Q,Q), g{CsXt) = /i(G,Ct) < /i(Cs,Ct) for (CO) ^ dU r\ I 
and thus there exist a sufficiently small neighborhood [/' C C of (Cs, Ct) and ((^^, Q[) G [/' such that 
0) < Ci') < fiiCs, Ct), which implies that d(s*, t*) = g{Q,Q) = 9{Cs, Ct) < oiCs, O- 

Hence, both /i and /2 are constant on C Then, by Lemma|7J ui,U2, and s* are coUinear and vi, 
V2, and t* are coUinear. Since m < 3, this situation violates the second part of Lemma|2] Thus, we get 
a contradiction again, concluding that m = |n(s* ,t*)\ > 3 for Case (BB) when both s* and t* lie on B. 

Case (BI): When s* £ B and t* e intV. We assume that s* e Cs e E and t* G intV. Define s(Cs) 
as done in Case (BB) with s{C*) = s*. We redefine the function /j : A — > M to be 

fi{Cs,tx,ty) := hl{s{Cs))+d{ulv'^ + hl{U,ty), 

where A := s-i((VP(it9 U VP(uj)) n e^) x (VP(v9 U VP(i;i)) is a subset of M^. Let g{Cs,t^,ty) ■= 
mini fi{Cs,tx,ty) for any {Cs, tx,ty) G flA- 

Analogously to Lemmas |3] and |5j each fi is convex on any convex subset of Di and there exists 
a convex neighborhood C C M'^ of {Q,t*,t*) with C C f]Di such that the four requirements are 
satisfied. Suppose that m = |n(s*,t*)| < 4. Then, Theorem [J implies that there exists a line ^ C 
through {Q ,t*,t*) G such that g is constant on £ n C, thus at least one of the fi is constant on ^ n C 
by Lemma m 

If only one of the fi is constant on ^ n C, then we have a contradiction as done in Cases (II) and 
(BB). Thus, assume that /i and /2 are constant on £ D C. By Lemmas |4] and |7J vi, V2, and t* should 
be coUinear. Further, by the second part of Lemma |2l t* must lie in the interior of the convex hull of 
Vt* . In order to have an interior point of the convex hull on the segment between vi and V2, we need at 
least two more points. Nonetheless, we have \Vt*\ < |n(s*,t*)| < 4, a contradiction. Thus, we have 
m = |n(s*, > 4 for Case (BI), as claimed. 

Finally, we complete a proof of Theorem [2l The claimed bounds on | Vg* \ and | Vf \ are shown by 
Lemma[2l and the bounds on |n(s*, are shown case by case as above. B] 

5 Computing the Geodesic Diameter 

Since a diametral pair is in fact maximal, it falls into one of the cases shown in Theorem |2l In order to 
find a diametral pair we examine all possible scenarios accordingly. 

Cases (V-), where at least one point is a comer in V, can be handled in 0{n^ log n) time by comput- 
ing SPM(i;) for every v V and traversing it to find the farthest point from v, as discussed in Section[2l 
We thus focus on Cases (BB), (BI), and (II), where a diametral pair consists of two non-corner points. 

From the computational point of view, the most difficult case corresponds to Case (II) of Theorem|2j 
in particular, |n(s*,t*)| = 5 in which 10 corners of V are involved, resulting in \Vs*\ = \Vt*\ = 5 
(see Appendix A.3). Note that we do not need to take a special care for the case of |n(s*,t*)| > 5. 
By Theorem |2] and its proof, it is guaranteed that there are five distinct pairs {ui,vi), . . . , (us, v^) of 
corners in V such that len^.^^. (s*, t*) = d{s*,t*) for any i G {1, • • • , 5} and the system of equations 
len„j^i,^(s, t) = • • • = lenMg,„5(s, t) indeed determines a 0-dimensional zero set, corresponding to a 
constant number of candidate pairs in intV x intT^. Moreover, each path-length function len„ ^, is an 
algebraic function of degree at most 4. Thus, given five distinct pairs (uj, Vi) of corners, we can compute 
all candidate pairs (s, t) in 0(1) time by solving the systemjl Then, for each candidate pair we compute 
the geodesic distance between the pair to check its validity. Since the geodesic distance between any 
two points s,t ^ V can be computed in 0{n log n) time [12], we obtain a brute-force 0{n^^ log n)-time 
algorithm, checking 0{n^^) candidate pairs obtained from all possible combinations of 10 comers in V. 

^Here, we assume that fundamental operations on a constant number of polynomials of constant degree with a constant 
number of variables can be performed in constant time. 
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As a different approach, one can exploit the SPM-equivalence decomposition of V, which subdivides 

V into regions such that the shortest path map of any two points in a common region are "topologically 
equivalent" fj). It is not difficult to see that if (s, t) is a pair of points that equalizes any five path-length 
functions, then both s and t appear as vertices of the decomposition. However, the currently best upper 
bound on the complexity of the SPM-equivalence decomposition is 0{n^'^) Q, and thus this approach 
hardly leads to a remarkable improvement. 

Instead, we do the following for Case (II) with iVg* \ = 5. We choose any five corners ui, . . . , U5 G 

V (as a candidate for the set Vg*) and overlay their shortest path maps SPM(iij). Since each SPM(Mj) 
has 0(n) complexity, the overlay consists of O(n^) cells. Then, any cell of the overlay is the intersection 
of five cells associated with wi, . . . , ^5 G F in SPM(ni), . . . , SPM(n5), respectively. Choosing a cell of 
the overlay, we get five (possibly, not distinct) w 1 , . . . , f 5 and thus a constant number of candidate pairs 
by solving the system len^^,^^ (s, t) = • • • = len^g^^g (s, t). We iterate this process for all possible tuples 
of five comers mi, . . . , us, obtaining a total of 0{v7) candidate pairs in 0{rJ log n) time. Note that the 
other subcases with \ Vs* \ < 4 can be handled similarly, resulting in 0{n^) candidate pairs. 

In order to test the validity of each candidate pair {s,t), we check the geodesic distance d(s,t) 
using a two-point query structure of Chiang and Mitchell [7 |: for a fixed parameter < 5 < 1 and any 
fixed e > 0, we can construct, in 0(n^+^*^^+'^) time, a data structure that supports ©(n^^*^ log n)-time 
two-point shortest path queries. Then, the total running time is 0(n^ log n) + 0{rv'^^^^~^'^) + 0{'n7) x 
0(n^~'^ logn). We set (5 = ^ to optimize the running time to 0(n^"'"n"'"'^). 

Also, we can use an alternative two-point query data structure whose performance is sensitive to the 
number h of holes Q: after Oin^) preprocessing time using O(n^) storage, two-point queries can be 
answered in 0(log n + h) timeO Using this alternative structure, the total running time of our algorithm 
becomes 0(n^(logn + h)). Note that this method outperforms the previous one when h = O(nTT). 

The other cases can be handled analogously with strictly better time bound. For Case (BI), we 
handle only the case of |n(s*,t*)| = 4 with |Vj*| = 3 or 4. For the subcase with |Vi*| = 4, we 
choose any four comers from y as , . . . , ^4 as a candidate for Vt* and overlay their shortest path maps 
SPM(fi). The overlay, together with V, decomposes dV into 0{n) intervals. Then, each such interval 
determines ui, . . . as above, and the side G on which s* should lie. Now, we have a system 
of four equations on four variables: three from the corresponding path- length functions len^^.t;, which 
should be equalized at {s* ,t*) and the fourth from the supporting line of e^. Solving the system, we get 
a constant number of candidate maximal pairs, again by Theorem |2] and its proof. In total, we obtain 
O(n^) candidate pairs. The other subcase with \ Vt*\ =3 can be handled similarly, resulting in O(n^) 
candidate pairs. As above, we can exploit two different structures for two-point queries. Consequently, 
we can handle Case (BI) in 0(n^^ 11 ^'^) or 0(n^(log n + h)) time. 

In Case (BB) when s*,t* G B, we handle the case of |n(s*,i*)| = 3 with \ Vs*\ = 2 or 3. For 
the subcase with \Vs*\ = 3, we choose three comers as a candidate of Vg* and take the overlay of 
their shortest path maps SPM(nj). It decomposes dV into 0(n) intervals. Then, each such interval 
determines three corners ^1,^2, f 3 forming Vt* and a side G on which t* should lie. Note that we 
have only three equations so far; two from the three path-length functions and the third from the line 
supporting to et- Since s* also should lie on a side Cg € E with eg 7^ et, we need to fix such a side Cg 
that PIkiO VP(nj) intersects eg. In the worst case, the number of such sides Sg is Q{n). Thus, we have 
O(n^) candidate pairs for Case (BB); again, the other subcase with \ Vg*\ = 2 contributes to a smaller 
number O(n^) of candidate pairs. Testing each candidate pair can be performed as above, resulting in 
0(n^+TT+'^) or 0(n^ (logn + h)) total running time. 

For Case (BB), however, one can exploit a two-point query structure only for boundary points 
on dV. The two-point query structure by Bae and Okamato [5] indeed builds an explicit representa- 

''if h is relatively small, one could use the structure of Guo, Maheshwari and Sack 1101 which answers a two-point query in 
O(ftlogn) time after 0{n^ ^ogn) preprocessing time using 0{n^) storage, or another structure by Chiang and Mitchell (7) 
that supports a two-point query in 0{h log n) time, spending 0{n + h^) preprocessing time and storage. 
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tion of the graph of the lower envelope of the path-length functions len„ restricted on dV x dV in 
0(n^ lognlog* n) time0 Since |n(s*,t*)| > 3 in Case (BB), such a pair appears as a vertex on the 
lower envelope. Hence, we are done by traversing all the vertices of the lower envelope. 
The following table summarizes the discussion so far. 



Case 


Independent of h 


Dependent on h 


(VV), (VB), (VI) 


0{'n? logn) 


(BB) 


0{rt' log n log* n) 


{ft' {log n + h)) 


(BI) 


0(n5+TT+^) 


0{n^{\ogn + h)) 


(11) 




0{n'^ {log n + h)) 



As Case (II) being a bottleneck, we conclude the following. 
Theorem 3 Given a polygonal domain having n comers and h holes, the geodesic diameter and a 
diametral pair can be computed in 0(n^"'"TT+'^) or 0(n^(log n + h)) time in the worst case, where e is 
any fixed positive number. 

We can avoid some difficult cases when /i is a small constant based on a simple observation: If there 
are two distinct shortest paths between s and t in V, then we know that there is at least one hole in the 
region closed by the two paths. In general, if k shortest paths exist between any two points of V, we can 
conclude /i > /c — 1. By contraposition, if h < k — 1, then there cannot exist two points that have k or 
more distinct shortest paths between them. 

Theorem 4 Given a polygonal domain having n comers and h holes, the geodesic diameter and a 
diametral pair can be computed in the following worst-case time bound, depending on h. 

• 0{n) time, ifh = (by Hershberger and Suri 071/ ). 

• 0{'n? logn) time, ifh= 1, 

• 0{n^ log n) time, ifh = 2 or 3, 

• 0(n''(logn + h)) time, if 4 < /i = 0{n^), 

• 0{'nJ^ri+'^'j time, otherwise. B] 



6 Concluding Remarks 

We have presented first algorithms that compute the geodesic diameter of a given polygonal domain. 
They are based on our new geometric observations on local maxima of the geodesic distance function, 
which show tight lower bounds on the number of shortest paths between a maximal pair. It is worth 
noting that with analysis in Section |5] the number of geodesic-maximal pairs is shown to be at most 
0{'nJ). On the other hand, one can easily construct a simple polygon in which the number of maximal 
pairs is O(n^). An interesting question would be what the maximum possible number of maximal pairs 
in a polygonal domain is. 

Though in this paper we have focused on exact geodesic diameters only, an efficient algorithm for 
finding an approximate geodesic diameter would be also interesting. Notice that any point s G "P and its 
farthest point t yield a 2-approximate diameter; that is, diam(T') < 2 maxtg-p d(s, t) for any s £ V. 
Also, based on a standard technique using a rectangular grid with a specified parameter < e < 1, one 
can easily obtain a (1 + e) -approximate diameter in 0((4 + ^) logn) timeO Breaking the quadratic 
bound in n for the (1 + e)-approximate diameter seems a challenge at this stage. Thus, we pose the 
following problem: For any or some < e < 1, is there any algorithm that finds a (1 + e)-approximate 
diametral pair in 0{n?~^ ■ poly(l/e)) time for some positive (5 > 0? 

''More precisely, in 0(n*A65(n) logn) time, where Am(n) stands for tlie maximum lengtli of a Davenport-Schinzel se- 
quence of order m on n symbols. 

^We thank Hee-Kap Ahn for giving us the idea of this approximation algorithm. 
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APPENDIX 



A More Examples and Remarks 

In this section, we show more constructions of polygonal domains and their diametral pairs with inter- 
esting remarks. In the figures, we keep the following rules: the boundary dV is depicted by dark gray 
segments and the interior of holes by light gray region. A diametral pair is given as {s* ,t*) and shortest 
paths between s* and t* are described as black dashed polygonal chains. 

A.l Examples where at least one point of a diametral pair lies on &P 




(d) (e) (f) (g) 

Figure 4: (a-c) Polygonal domains whose geodesic diameter is determined by a comer s* and (d-g) variations 
of the construction (c). (a) When both s* and t* are corners; (b) When t* is a point on &P; (c) When t* e iniV. 
This polygonal domain consists of two holes, forming a narrow corridor and three shortest paths between s* and 
t* . Here, we have d(s*, wi) = d{s* ,V2) = d(s*, W3) and t* is indeed the vertex of SPM(s*) defined by Vi,V2,v^\ 
(d) Variation of (c) with all convex holes; (e) Three shortest paths are not enough to determine a boundary- 
interior diametral pair; (f) If we add one more hole, then the diameter is determined by s* e S and t* e intT^ 
with four shortest paths; (g) A polygonal domain made by attaching two copies of (e) and modifying it to have 
d(wi, vi) = d(M2, i'2) = d(M3, U3). Observe that, in this polygonal domain, the diameter is determined by two 
boundary points with three shortest paths. 

Note that, as expected, every example in Figure|4]satisfies Theorem|2l An interesting construction is 
FigureHtg), where neither of the two centers of Anin2'U3 and of Ai;if 2^3 appears in any diametral pair. 
Also note that Figure Ud) consists of convex holes only. We think that any complicated construction 
can be "convexified" in a similar fashion. This would suggest that computing the diameter in polygonal 
domains with convex holes only might be as difficult as the general case. 
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A.2 A proof for Figure [He): Case (II) with 6 shortest paths 




Figure 5: A schematic diagram corresponding to the polygonal domain shown in FigurefTJc). 

Claim 1 In the polygonal domain described in FigureW[c), {s*,t*) is the unique diametral pair. 
Proof of Claim. Recall that by construction of the problem instance, the triangles AU1U2U3 and AviV2V^ 
are regular and d{ui,vi) = d{ui,V2) = d(n2, ^2) = d(n2, wa) = d(n3, vs) = d{us, vi) = L, for some 
arbitrarily large value L > 0. Also, s* and t* are the centers of AuiU2Us and AviV2Vs, respectively. 

We assume that both triangles AuiU2Us and AV1V2V3 are inscribed in a unit circle. Then, we have 
d{s*,t*) = 2 + L. For any point s on any shortest path between Ui and Vj, it is easy to see that 
d(s, t) < -v/3 + L < d{s*,t*) for every point t £ V. Thus, any point on those paths cannot contribute 
to the diameter. 

(1) First, observe that max(gAt)ii;2f3 d(s*, t) = maXsgAMiMa^s ^i^? **) = d(s*, t*). 

(2) For any s G AU1U2U3, its farthest point t G AviV2Vs is on the angle bisector of some Vi. 
Consider any s G AU1U2U3. Without loss of generality we assume that ||s — < mini{||s — Ui\\}. 
Then, both the shortest paths to vi and to V2 from s pass through ui. We thus have d(s, vi) = d(s, ^2) 
by construction and its farthest point t G AV1V2V3 must be in the angle bisector of f 3. By symmetry, the 
same property holds when the closest corner from s* is either U2 or ^3. 

Conversely, for any t, its farthest point s G AU1U2U3 must be on a bisector of some n,. In any 
diametral pair (s, t), we have that t is the farthest point of s (and vice versa), so both must be on one of 
the angle bisectors. 

(3) If (s, t) is a diametral pair, then s G Ui-s* and t G Vjt*, for some i and j. Suppose that s lies on 
the bisector of ui but not in between ui and s*. We then have \\s — U2\\ = \\s — u^W < \\s — ui\\ and 
d(s, vi) = d{s, V2) = d(s, ^3) = ||s — U2\\+Lby construction. This implies that t* is the farthest point 
of such s. Since ||s — M2II < 1 and thus d(s, t*) < 2 + L, (s, t*) is not a diametral pair. 

(4) Now, pick any point s G uis* with s 7^ s*. Suppose that t G AV1V2V3 is the farthest point from 
s. We know that t G v^t* by above discussions. In this case, we have four shortest paths between s 
and t through (ui, vi), (ni, ^2), (m2, t's), and (^3, ^3); the other two are strictly longer unless s = s*. 
Thus, by Theorem|2j such s G uis* with s 7^ s* and its farthest point t cannot form a maximal pair. By 
symmetry, the other cases where s £ UiS* can be handled. 

Hence, (s*, t*) is a unique diametral pair and the geodesic diameter is 2 + L. 
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A.3 Diametral pair of Case (II) with exactly 5 shortest paths 

Here, we present a polygonal domain in which the diameter is determined by two interior points and 
exactly five shortest paths between them. This proves the tightness of Case (II) in Theorem |2l 




Figure 6: A schematic diagram of a polygonal domain in which | V,. | = \Vt* \ = 3 and |n(s*, = 5. 

Figure |6] shows a schematic description of a polygonal domain V. We assume that only the position 
of the vertices Ui and the Vi are geometrically precise. Thus, we construct the problem instance such that 
we have ui = U2 = U3, vi = v^, and = v^, and the convex hulls of the Uj and of the Vi form isosceles 
triangles A„ and A„. Each of and A^, is inscribed in a unit circle centered at c„ and Ou. Moreover, the 
bases of both triangles are horizontal and the angles opposite to the bases are 18° and 112°, respectively. 
Note that the side lengths of the triangles A^ and A^ are as follows: \\ui — u^\\ = 1.97537- • • and 
||t(4 -U5II = 0.61803- - •; \\v2 - vi\\ = 1.11833 - •• and \\vi - V3\\ = 1.85436- • •. 

In this configuration, we set the constants as follows: letting L := d{ui,vi) = d(u3, v^) be some 
sufficiently large number, we set d{u2, V2) = L + 0.5 and d(u4, ^4) = d(u5, v^) = L + 0.2. Note that 
this configuration can be realized with four obstacles in a similar way as Figure [He). 

Since we have fixed all necessary parameters, we have a fully explicit description of the leiiu^ Due 
to the difficulty of finding an exact analytical solution, we used numerical methods to solve the system 
of equations len^^ „j (s, t) = • • • = leuu^^y^ (s, t). We have found that there is a unique solution {s*,t*) 
such that s* G A^andt* G A^ ; we obtained s* = c„ + (0, -0.102795 •••), = c„ + (0, 0.555361 ••• ) 
and d{s*,t*) = 2.047433734 ■ ■ ■ + L. (See Figure[6l) 

We first checked that {s*,t*) is indeed a maximal pair based on the following lemma, which can be 
shown using elementary linear algebra together with the convexity of the path-length functions. 

Lemma 8 Suppose that (s, t) is a solution to the system lenu^^vi {s,t) = ■ ■ ■ = len^g^^^ (s, t). If any four 
of the five gradients Vleiiu^^^. at {s, t) are linearly independent (as vectors in a A-dimensional space) and 
one of them is represented as a linear combination of the other four with all "negative" coefficients, then 
{s, t) is a local maximum of the pointwise minimum of the five functions len„. i,. . H 

Next, to see that {s*,t*) is a diametral pair, we have run our algorithm for each of Cases (BB), (BI), 
and (II); as a result, there are 44 candidate pairs, including {s*,t*), falling into those cases among which 
at most 11 are maximal and only {s*,t*) is diametral. Note that the pair (s*, t*) is the only candidate 
pair of Case (II). Also, observe that any point on the shortest path between m and Vi cannot belong to a 
diametral pair. This implies that none of the Uj and the Vi belongs to a diametral pair and thus that there 
is no diametral pair in Cases (V-). In addition, we also sampled about 350,000 points uniformly from 
each of A^ and A^,, and evaluated the geodesic distances of the 350,000^ pairs. 

Note that one can modify the construction to have \ Vs* | = | Vt* | = |n(s*, = 5. For the purpose, 
we can split ni, ^2, U3 into three close comers (analogously for corners, ^1,^5 and V3, V4). The splitting 
process should preserve the differences between the distances d{ui,Vi) for alH = 1, . . . , 5. We also 
have tested such an example in the same way as above and concluded that a solution equalizing the five 
path-length functions is indeed a diametral pair. 
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B Proofs Omitted from the Body 



Proof of Lemma [U We only give a proof for the second statement, which imphes the first. The case of 
d = lis trivial, so we assume d > 1. Let Hi, . . . , Hm be any m closed hemispheres on S^^^, and hi be 
the hyperplane through the origin in such that Hi lies in a closed half-space supported by hi. In this 
proof, we denote by Hi the open hemisphere, defined to be Hi = Hi \ hi. Also, let T-Lj := ni<j<j 
andnj := r\i<i<j Hi. 

Suppose that Tim = 0- Let k be the smallest integer such that Tik = O- By definition, k > 2 and 
'Hk-i 7^ 0- Since the intersection of any k — 1 non-parallel hyperplanes of includes a (d — A; + l)-flat 
and each hi contains the origin, ni<i<fc-i ^« includes a. {d — k + l)-flat through the origin and thus 
Tik-i includes a great {d — A;)-sphere G on S'^^^. Since x £ G implies —x G G for any x G S'^^^, we 
must have G C h^, in order to have an empty intersection Tik- This implies that ni<i<fc ^Iso includes 
a. {d — k + l)-flat through the origin, and further that ni<i<m includes a (d — m + l)-flat through the 
origin. We hence conclude that T-Lm = ni<j<m ^i includes a great {d — m)-sphere on S"^^^. 13 



Proof of LemmalU Since {s*,t*) is a maximal pair, ds* (t) := d{s*,t) is maximized at t* over a small 
neighborhood U C V of t*. Thus, as discussed in Section|2l if t* ^ V, then t* must be either a vertex of 
SPM(s*) or an intersection point between an edge of SPM(s*) and dV. If t* e intV, then t* should fall 
into the former case and we have at least three corners vi,V2,V3 G V determining a vertex of SPM(s*). 
If t* £ B, then t* may fall into either case. Even in the latter case, t* must lie on B and we have at least 
two corners 1^1,^2 G V determining an edge of SPM(s*). 

The second part of the lemma can be shown as follows. If t* S intT^ but lies out of the interior of 
the convex hull of Vt* , then we can find another point t £ U for any neighborhood U C infP of t* 
such that \\t - fill > \\t* - Vi\\ for every Vi G Vf, implying that d(s*,i) > d(s*,t*). \ft*£e£E 
but every Vi G Vf lies on the supporting line i of e, then when we move t* in a perpendicular direction 
to I, we obtain a strictly larger distance as above. (Remark that a similar argument can be found also 
in HSl Lemma 2.2].) S 



Proof of Lemma |3] Since the sum of convex functions is a convex function, it suffices to show that hf 
and h\ are convex. More precisely, for any (si, ti), (§2, ^2) G C and < A < 1, we have 

/,(A(si, ti) + (1 - \){S2M)) = KiXsi + (1 - A)s2) + dK, v',) + hiiXti + (1 - X)t2) 

< xh^isi) + (1 - x)ht{s2) + dK, + xhlih) + (1 - x)hi{t2) 

= Xfi{si,h) + il-X)f,{s2,t2) 

if hf and hj are convex. 

We thus show the convexity of hf on any convex subset Cs C VP(ti^) U VP(ui). The convexity of 
hj can be shown in the same way. There are two cases: u'^ = Ui or u[ ^ Ui. For the former case, /i| is 
convex on Cg since it measures the Euclidean distance between Ui and a given point in Cg. Letting be 
the line through Ui, u'^, and also s*, for the latter case, Cg may be partitioned by £q into two regions Ai 
and A2, where Ai = GgC] VP(n^) and A2 = Cs \ Ai. Note that /if is convex on Ai and on A2. Thus, 
we are done by checking every point on ^0 H C^. 

Pick any s G £0 H Cs and any line ^ C through s. Let 6 be the angle between and I. If we 
restrict the domain of /if on £ n Gg, then one can check with basic calculus that both the derivatives of 
||s — Will + ||nj — n^ll and of ||s — ii^|| are equal to ccosO at s for some constant c. Hence, hf is smooth 
and convex along I. Since we have taken any line I through any point on £9 H C^, this suffices to prove 
the convexity of /if on Cs. B] 
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Proof of LemmalU There are four cases according to the position of s and t: either s G VP(u^) or 
s G VP(ni) \ VP('u.); either t £ VP(v-) or t G \/P{vi) \ \/P{v[). We give a proof only for the case of 
s G VP(tii) \ VP(m-) and t G \/P{vi) \ VP(t>-) but the proofs for the other cases are almost identical. 

Recall that s Ui and t ^ Vi. Any ray (or half-line) 7 C M"^ with endpoint (s, t) can be determined 
by three parameters {Og, 6t, A) with < 6s,0t < tt and A > as follows: Let 7^ and 7^ be the projections 
of 7 onto the [s^, Sy)-plane and the {tx,ty) -plane, respectively. Note that 7^ is a ray in the {sx, Sy)-plane 
with endpoint s and jt is a ray in the {tx,ty) -plane with endpoint t. Let Og be the smaller angle at s 
made by 7^ and another ray starting from s in direction away from Ui. Define 9t analogously. Then, 
the derivative of /j at (s, t) along 7 is represented as c(cos 6*^ + A cos 9t) for some constants A > and 
c > 0. Also, the second derivative of at (s, t) along 7 is derived as c (pzj^ + Ap^^) . 

Suppose that fi is constant along 7 locally around (s, t). Then, its first and second derivatives should 
be zero in a small neig hborhood [/ C of {s, t) with U C Di. First, we observe that A should be 
positive; if A = 0, then t is fixed while s moves along 7^ so that fi does not stay constant. Since every 
term of the second derivative is nonnegative and A > 0, we get {Os^Ot) = (0, vr) or (vr, 0). We hence 
obtain only two solutions {Os, 6t, A) = (0, vr, 1) or (vr, 0, 1). The two rays 7 corresponding to these two 
solutions form a unique line £i C through (s, t) such that fi is constant on £i n U. 

Now, we pick any (s', t') G £j fl intDj such that both (s, t) and (s', t') lie on a common connected 
component of ii n intDj. Then, we observe that fi{s', t') = fi{s, t). Since any convex neighborhood C 
of (s, t) with C <Z Di intersects li in a single connected component, for any (s', t') G n C, it holds 
that fi{s',t') = fi{s,t), proving the second part of the lemma. 

Remark that £i n intDi may consist of more than one connected components; a typical situation 
where ii n intDi is disconnected is when £i passes through {ui,t") or {s", Vi) on dDi for some s" G 
\/P{u'j) U VP(uj) and some G VP(i'j') U VP(i;j). See Figure |3] for more intuitive and geometric 
description on ^j. H 



Proof of LemmalU (i) In this proof, we define len„^„(s,t) = 00 if s ^ VP(u) or t ^ VP(f). By 
definition, there exists a small neighborhood Ui C Di of {s*,t*) such that fi{s, t) = min{len„-^^- (s, t), 
len„^_„.(s,t),len„.^^^(s,t),lenu/^„/(s,t)} = minu^^g^.ny len„,t,(s, i) for all G Ui. We claim that 
there exists a convex neighborhood C C f]iUi such that for any (s, i) G C 

d(s,i)= mill fi{s,t) = g{s,t). 

l<t<m 

To prove our claim, assume to the contrary that for every convex neighborhood C C of {s*,t*) G 
there exist a pair (u, t>) of corners and (s, t) £ C such that d(s, t) = lenu,v{s, t) < mirij /j(s, t). 

Note that none of the shortest paths vTj G n(s*, t*) between s* and i* passes through both of such u and 

V since otherwise lenu^^(s, t) = minj t). 

Consider a sequence Ci, C2, ... of neighborhoods of (s*, t*) G that converges to the singleton 

{(s*, t*)}. Since there are only -n? pairs of corners, there exist a fixed pair (uq, fo) of corners and a 

subsequence C^^, C^j, • • • converging to the singleton {(s*, t*)} such that none of the vTj passes through 

both of no, vq and for any integer j > there is a point {sj,tj) G Ck^ with 

d{sj,tj) = lenuo,vo{sj,tj) < min fi{sj,tj). 

l<i<m 

Since limj_>,oo('Sj, tj) = {s* ,t*), it holds that linij^oo d{sj,tj) = limj^oo minj /i(sj, tj) = d{s* ,t*). 
By the sandwich theorem, we thus have 

lim lenuo,vo{sj,tj) = len„ ^ (s*,t*) = d{s*,t*). 

This implies the existence of the (m+l)-st shortest path between s* and t* since none of the tTj G 
n(s*, t*) contains both uq and vq, a contradiction. 
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(ii) The claim follows from the fact that fi{s*, t*) = len„,,„, {s*,t*). 

(iii) From (i), g{s, t) = d{s, t) for all (s, t) € C. Since d attains a local maximum at (s*, t*), so does g. 

(iv) Consider any pair {s,t) such that s G s*Ui and t € but s 7^ s* and t 7^ t*. Then, fi{s,t) = 
len„._„.(s,t) < fi{s*,t*). Hence, there is no neighborhood [/ C of {s*,t*) such that fi{s*,t*) < 
fi{s, t) for any (s, t) € ?7, implying that {s* ,t*) is not a local minimum of any /j. H 



Proof of Lemma|6] Let li C be the unique line through {s* ,t*) such that fi is constant on ii n C. 
The uniqueness of is proven by Lemma |4] 

Observe from the proof of Lemma |4] that the projection of £i onto the (s^, , Sy)-plane is the line 
through s* and Ui. Also, the projection of li onto the {tx, ty)-plane is the line through t* and Vi. Hence, 
li = ij implies that Ui, uj, s* are collinear and Vi, Vj, t* are coUinear. First, since the pairs (uj, Vi) are 
all distinct, we have Ui 7^ Uj or Vi ^ vj. If Uj = Uj and Vi 7^ f j, one can easily check that £i 7^ £j. We 
thus have Ui 7^ Uj and Vi 7^ Uj. Moreover, s* must lie in between Ui and Uj and must lie in between Vi 
and Vj by definition; if lies in between Ui and s*, then the first comer of vTj from s* becomes Uj since 
the three are collinear. Therefore, for each i G {!,..., m}, there is at most one index j G {1, . . . , m} 
such that i j and = £ j , completing the proof. B] 



Proof of Lemma |7] The lemma immediately follows from Lemma |4] and its proof. B] 
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